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Abstract 

This article gives a contemporary and to some extent pedagogical review of the 
current theoretical understanding of the formation of the superconducting state in 
metallic systems with a variable density of carriers. We make an attempt to de- 
scribe the crossover from the Bose-Einstein condensation type (small densities) to 
the Bardeen-Cooper-Schrieffer one (large densities). The functional methods are 
used throughout the treatment. The most of the results are considered in a review 
form for the first time. Some of them (in particular, the possible opening of a pseu- 
dogap) are used to explain the experimental data available for high-temperatures 
superconductors. 
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1 Introduction 



1.1 The general formulation of the problem 

It is well known that there are two fundamental limiting descriptions which allow us to un- 
derstand the equilibrium low temperature " superphenomena" , such as superconductivity 
and superfluidity, in 3D systems: 

i) Bardeen-Cooper-Schrieffer (BCS) theory [l], in which the normal state is a degenerate 
Fermi liquid that undergoes a cooperative pairing (Cooper) instability at the temper- 
ature Tp. Here two processes, the formation of Cooper pairs and their condensation 
(macroscopic occupation of a single quantum state), occur simultaneously at the transi- 
tion temperature, T^ D = T P . 

ii) Bose-Einstein condensation (BEC) of bosons (see e.g. 0) in another single quantum 
state at the temperature Tb- In fact, these bosons (for instance, 4 He) are tightly bound 
(composite) particles made up of an even number of fermions. The bosonic particles are 
formed at some high temperature Tp of their binding (dissociation). However in contrast 
to the BCS case, these "pre-formed" bosons condense only at T B <C T P . The above- 
mentioned single quantum states as well as an intermediate case, which we will discuss 
below in details, are usually described by an order parameter which is homogeneous in 
the 3D case assuming that there are no external fields. 

Most, if not all, of the real 3D systems clearly fall into either category i) or ii). 
For example, 3 He and essentially all metallic superconductors that we now understand 
are Fermi superfluids described by the first category, whereas 4 He is a Bose superfluid 
described by the second one. 

Nevertheless, the general problem of the crossover (or interpolation) from the BCS 
scenario of superconductivity with cooperative Cooper pairing to the formation of com- 
posite (separate) bosons and their BEC has remained of great interest for a long time. It is 
connected primarily with a deeper understanding of the phenomenon of superconductivity 
even for 3D systems which are clearly closer to the BCS limit. 

The significance of the problem increases when one tries to reduce the dimensionality 
of space from 3D to 2D. In such a case the problem acquires not only fundamental sig- 
nificance, but also practical importance. The latter, as we shall discuss below, is mainly 
related to the discovery of the high-temperature superconductors (HTSC) ||. It is not 
the only reason for the widespread interest in the problem of the crossover. Even from 
the theoretical point of view the 2D crossover is far more mysterious than the 3D one. At 
first glance the BCS limit with weak attraction and high fermion density is completely 
different from the Bose limit with strong coupling and low boson density. Indeed, the 
simplest BCS theory predicts a finite value of T C 2D = TjP while the BEC is forbidden in 
2D systems for massive bosons @ for which = 0. This is the reason one might ask 
whether it is possible to obtain a 2D crossover. 

It is worthwhile to note here that by virtue of the Coleman- Mermin-Wagner-Hohenberg 
theorem M a "super-behaviour" [] with homogeneous order parameter (the latter related 
to breaking of some continuous symmetry) is not possible for pure 2D systems. This is 

Super-behaviour" refers to 2D ferro- or antiferromagnetism in the case of the Mermin- Wagner the- 
orem and to superfluid behaviour in Hohenberg's theorem, while the Coleman theorem is the field- 
theoretical generalization of the previous theorems (see e.g. [Q). 
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due to the fluctuations of the order parameter (in particular of its phase) destroying the 
long-range order (see also 0). Therefore, the temperature T C 2D of establishing long- 
range order has to be zero. The statement that T| D = is in agreement with this general 
theorem. In fact the BCS result Tp D does not contradict these theorems if one remembers 
that as a rule the BCS theory presupposes the mean field (MF) approach. Thus it does 
not take into account any fluctuations and the result only means that T C MF ^ = T 2 P D + 0, 
while the temperature T C 2Z? of actually establishing the long-range order is zero in this case 
as well. Such a conclusion, however, results in a general question: how can one study 2D 
models of superconductors with finite T c values? 

There are at least two ways to answer this question. 
First, 2D models are only a mathematical idealization and any real system must be at least 
quasi-two dimensional (quasi-2D). This immediately gives T c ^ 0, although T c < T^ F . 
Second, there is another interesting possibility for 2D systems to pass into the superfluid 
state at a temperature which is usually denoted by Tbkt- This corresponds to the so called 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition. The BKT superfluid state has an 
inhomogeneous order parameter and is therefore not forbidden by the above-mentioned 2D 
theorems. The formation of this state is quite different from the ordinary superconducting 
transition and results in a new way of condensation (see Chapter 5). 

The possibility for real systems to undergo the BKT transition crucially depends on 
their spatial anisotropy. Thus, one should study both 2D and quasi-2D cases. 

Another important aspect of the crossover problem is related to the following unusual 
property of the HTSC compounds. Almost all the physical characteristics of HTSC's (T c 
included) crucially depend on the itinerant carrier density. This means that the crossover 
in superconducting systems cannot be studied without taking into account the changing 
delocalized fermion density. We believe that many of the HTSC "anomalies" are caused by 
the unusual (two-stage) formation of the superconducting state in 2D and (very probably) 
quasi-2D systems with an arbitrary carrier density 

Hopefully, this short introduction has convinced the potential readers that the su- 
perconducting condensate formation at the different densities of interacting fermions (in 
particular, in the BCS and Bose limiting cases) deserves a separate review. 



1.2 History 

The idea that the composite bosons (or, as they called, local pairs) exist and define the 
superconducting properties of metals is in fact more then 10 years older than the BCS 
theory. As early as 1946, a sensational communication appeared saying that the chemist- 
experimentator Ogg had observed superconductivity in the solution of Na in NH 3 at UK 
||10|| . It is very interesting that the researcher made an attempt to interpret his own result 
in terms of the BEC of paired electrons. Unfortunately, the discovery was not confirmed 
and both it and his theoretical concept, were soon completely forgotten (for the details 



see e.g. [p"T|j). 

It is appropriate to note that the history of superconductivity has many similar exam- 
ples. Probably, this explains why Bednorz and Miiller named their first paper "Possible 
high-T c superconductivity..." ||. And even today there appear many unconfirmed com- 
munications about room-temperature superconducting transitions. 

A new step, the development of the local pair concept, was taken in 1954 by Schafroth 
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who in fact re-discovered the idea of the electronic quasi- molecules [12]. This idea was 
further developed in the Schafroth, Blatt and Butler theory of quasi- chemical equilibrium 
li3f , where superconductivity was considered versus the BEC. Such a scenario, unfortu- 
nately, could not compete with the BCS one due to some mathematical difficulties which 
did not allow the authors to obtain the famous BCS results. Then the triumph of the 
BCS theory replaced the far more obvious concept of the local pairs and their BEC by the 
Cooper ones and their instability which takes place in the necessary presence of a Fermi 
surface (more precisely a finite density of states at the Fermi level). In contrast for the 
local (i.e. separated) pairs the formation is not in principle connected to this density of 
electronic states. Also unlike the local pairs, the Cooper ones are highly overlapping in 
real space. Therefore, the Cooper pairing should be understood as a momentum space 
pairing. 

Later experiment reminded us about the local pairs. It was Frederikse et al. 0] who 
found that the superconducting compound SrTiOs (T c ~ 0.3K) has a relatively low density 
of carriers and, moreover, this density is controled by Zr doping. The first deep discussion 
of the possibility of the BCS-Bose crossover and pairing above the superconductivity 
transition temperature for a low density of carriers was carried out by Eagles QX5j| (see 
also his et al. relatively recent article [TIJ) in the context of SrTiOs. 

The further history of the investigation of the BCS-Bose crossover has been often cited 
in the current literature. So, we note only that the features of 3D crossover at T = were 
considered in [ H]] , and its extension to finite temperatures was first given by Nozieres and 
Schmitt-Rink [18 |. 

The 2D crossover in superfluid 3 He was studied in 



IS], where the very natural and 



convenient physical parameter Eb, the bound pair state energy, was first used. 

The discovery of HTSC significantly stimulated the interest in the problem of crossover 
and related phenomena. There are many articles dealing with its study and we can only 

120, 



2T1, p2l pi E3, 



3 Ell 



sec 



mention some of them at this point 
also the recent excellent review devoted to other aspects of the BCS-Bose crossover, 
e.g. lattice models, time-dependent Ginzburg-Landau theory, numerical study of normal 
state crossover etc., which we because of this shall not discuss). We note only that some 
historical debates are associated with the parallel (and generally speaking independent) 
development of the macro- and microphysics of superfluidity and superconductivity since 
their discoveries were touched upon by Ginzburg in his new extremely comprehensive 
survey [B3|. 



1.3 Relevance to the high-temperature superconductors 

As we have already mentioned the crossover problem appears to be relevant to the general 
problem of the understanding of the HTSC. Indeed, these superconducting compounds 
have some peculiarities which place them much closer to the Bose or at least to the 
crossover region than majority of low temperature superconductors. We shall discuss the 
peculiarities in detail here. 

Certainly, it is necessary to point out right away that the crossover phenomena do not 
(and cannnot) address the problem of the mechanism for HTSC. The problem of HTSC 
itself is very difficult and controversial. Nevertheless the treatment of the crossover may 
shed light on some of the features of HTSC. 
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Almost all HTSC reveal following common properties: 



• a relatively low and easily varied density of itinerant (appeared owing to the doping) 
carriers; 

• a 2D, or more exactly quasi-2D, character of the conductivity and the magnetism; 

• the block structure of the crystal lattices, i.e. from 1 to 6 superconducting CuC>2 
layers per unit cell. 

Let us consider these items. The ground state of the copper-oxide based materials 
forms due to the strong antiferromagnetic spin fluctuations in the proximity of the metal- 
insulator transition (see e.g. review [53[). It seems plausible that further doping of these 



materials results in the appearance of weakly interacting (i.e. non-strongly correlated) 
itinerant carriers (holes). However the exact nature of the ground state including the 
strong electron-electron correlations is in fact not yet established. The density rif of these 
holes is not as large as in ordinary metals, so the mean distance between them proves 
to be comparable with a pair size £o or a coherence length 0. This situation is signifi- 
cantly different from the conventional BCS theory where the parameter £ greatly exceeds 

1 /2 

the mean distance between carriers which is ~ . Experimentally the dimensionless 
value of kp^Q (kp is the Fermi momentum) which describes the ratio of the pair size and 
the distance between carriers is about 5-20 for HTSC while for the low-temperature 
superconductors it is about 10 3 - 10 4 |20, |34 ]. 



It follows from the above that the new materials are likely to be in an intermediate 
regime between the Cooper pairs and the composite bosons, at least when the doping is 
not large and the value of T c is far from the highest possible (optimal) one. Although 



there is some evidence [Z^] that HTSC at the optimal doping are closer to the BCS limit 
than to the Bose one. 

Besides, because the coherence length is less than a lattice spacing in the direction 
perpendicular to the Cu0 2 planes (c-direction), the superconductivity in the copper ox- 
ides takes place mainly in the isolated Cu02 layers (or their blocks). This is the reason 
why pure 2D models of HTSC are commonly accepted. Undoubtedly these cuprate layers 
are connected even if they are situated in different unit cells. The mechanisms for this 
connection (by coherent or incoherent electronic transport) are not yet known exactly. 
There is no question, however, that one has to take into account the possibility of differ- 
ent (for instance, direct or indirect) interlayer hoppings to develop a full theory of HTSC. 
Therefore, strictly speaking, we need to consider quasi-2D models of these superconduc- 
tors. Moreover, as we have already mentioned, it is necessary to extend the system into a 
third- dimension in order to have a finite value for T c . It can be seen from the anisotropy of 
the conductivity that the influence of the third-dimension varies strongly from one family 
of cuprates to the other. For example, the anisotropy reaches 10 5 for the Bi- and Tl-based 
cuprates, while its value for YBaCus06+5 compound is close to 10 2 . It seems plausible 
that in the Bi(Tl)-compounds where the transport in the c-direction is incoherent that 
the scenario with BKT transition is more probable. In the Y-ones this transport is rather 



2 It is well to bear in mind that the coherence length strictly speaking is distinguished from the pair 
size, especially at low carrier density (this question will be treated in Section 3.4). 
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coherent and the superconducting state, at least at high doping, has a 3D character, with 
a homogeneous condensate which appears in the ordinary way. 

In general, it is now evident that, depending on the crystal anisotropies (and/or the 
relationship between the intensity of c-transport and the doping value) the quasi-2D 
systems can undergo two or even three phase transitions. The first one is from the normal 
phase to another normal (with a pseudogap) phase with a finite density of incoherent pairs. 
The second one is from the latter phase to the BKT phase with an algebraic order. The 
third and final one is from the BKT phase to the superconducting phase with long-range 
order. These possibilities will be discussed below. 

The role of interplanar effects inevitably increases when the superconducting planes 
are situated in the same unit cell so that the carriers in the adjacent layers are strongly 
coupled and can form interplanar pairs. This possibility together with the concentration 
effects brings about the simultaneous existence of the Cooper and local pairs in the system 
@- 

1.4 Outline 

The remainder of this review is organized as follows. In Chapter 2 we first describe the 
model for studying the BCS-Bose crossover, especially in 3D systems. Then in Section 2.2 
we introduce in detail the functional integral formalism which will be used throughout. 
The rest of this Chapter is devoted to applying the methods developed to the 3D crossover. 

The zero temperature 2D crossover is discussed in Chapter 3. Considerable attention 
will be given to the effective potential (see also Appendix A), which allows one to obtain 
the main equations. The kinetic terms of the effective action are obtained in Section 3.4. 
This gives one the possibility to investigate the dependencies of the coherence length and 
the penetration depth on the carrier density. 

In Chapter 4 the quasi-2D crossover is studied. 

Finally, the finite temperature 2D crossover is addressed in Chapter 5. We describe 
how it proves necessary to modify the accepted methods to consider the BKT phase 
formation (refer also to Appendices B and C) and its dependence on carrier density 
correctly. An attempt to explain some observable normal state anomalies HTSC is also 
made. 
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2 3D crossover: critical temperature 



It is most convenient for studying the crossover to start with the 3D case where there is no 
problem with establishing the long-range order below T^ D . We shall mainly follow here 
the paper |2(| (see also the review l32fl) in which the corresponding results were obtained 
by applying the more appropriate in this case the functional integral formalism rather 
than to the original paper 



1* 



2.1 Model 

Let us introduce a continuum field model of fermions with an attractive two-body inter- 
action. Our goal is to consider how the critical temperature T c 3£> chang function 
of attraction and to establish the main factors which determine it in the BCS and Bose 
regimes. 

The simplest model is described by the Hamiltonian density 

H = -^l{x) + f?j ~ ^}(x)Vi(x)^(x)V T (x), (2.1) 

where x = r, r (r is a 3D vector); ip a {x) is the Fermi field; m is the fermion effective mass; 
a =|, I is the fermion spin; V > is the attraction constant. The chemical potential 
fx fixes the average density nj of the free (bare) carriers. We choose units in which 
h = kg = 1 and the system occupies the volume v. 

Since we want to describe the Bose regime in terms of its constituent fermions, we 
have to allow the magnitude of the attraction V to be arbitrary. In addition, we cannot 
use a simple BCS-like cutoff uj d <C ep (usually wo is the Debye frequency) since we must 
allow the possibility for all the fermions to be affected by the interaction, not just for a 
small fraction (~ ud/cf) in a shell around the Fermi energy e^. We shall thus study a 
dilute Fermi gas in which the range of the attractive interactions can be characterized by 
a shape independent parameter, the scattering length a s in three dimensions. We shall 
describe in Section 2.3 how the " renormalized" coupling a~ x replaces "bare" V. 



2.2 Formalism 

As was already noted, the functional integral approach along with the Matsubara thermal 
technique is more appropriate for the problem that is studied here and in the subsequent 
Chapters. Thus, let us consider the formalism used. 

First of all, introducing Nambu spinors for fermion fields [^] (see also [[I]]) 

= * t(x) = ( WW ^(*)) ( 2 - 2 ) 

one has to rewrite ( |2.1| ) in a more suitable form: 

U = -^Ux) (— + u] T 3 ^(x) - V^(x)t + ^(xW(x)t^(x), (2.3) 
\2m J 

where r 3 , t± = (ti ±ir 2 )/2 are Pauli matrices. 
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Now the partition function is expressed through the Hamiltonian ( |2.3|) as: 



Z(v, fj,, T) = J V^VW exp l-J^drJ dr[^(x)d T ^{x) + H{r)\ 



(2.4) 



where (3 = 1/T and V^'D^^ denotes the measure of the integration over the Grassman 
variables \& and satisfying the antiperiodic boundary conditions: ^(r, r) = — ^(r + 
P, r) and *t(r, r) = -*t(r + /3, r). 

If it was possible to calculate the partition function ( |2.4| ) one could, in principle, obtain 
all thermodynamical functions from the thermodynamical potential 



Q(v,fi,T) = —T\nZ(v,fM,T). 



(2.5) 



Using now the auxiliary Hubbard-Stratonovich complex scalar field in the usual way 
one can represent G2.4|) in an equivalent form: 



Z{v,n,T) = J V^V^VQV®* exp |- £ dr J dr 
'V 2 



\$(x) 



d T I - r 3 



2m 



_$(xl 



V 



(2.6) 



The main virtue of this representation is a nonperturbative introduction of the composite 
fields $(x) = V&(x)t + V(x) = V^(x)$\(x), $*(x) = 

V^(x)t-^(x) = Vtpi(x)ip^(x) and the possibility to develop a consistent approach. 
Specifically, the expression ( |2.6|) turns out to be rather convenient for studying such 
a nonperturbative phenomenon as, for example, superconductivity. In this case the com- 
plex Hubbard-Stratonovich field naturally describes the order parameter arising due to 
the formation of the Cooper or the local pairs. The average value of |$|(= A) is propor- 
tional to the density of pairs, on the one hand, and determines the gap in the one-particle 
Fermi-spectrum, on the other. 

The integration over fermion fields in ( |2.6j ) can be done formally even though $ and 
$* depend on the spatial and temporal coordinates. Thus, one obtains (formally exactly) 



where 



Z(v,fi,T) = J P$P$*exp[-/30(u,/i,T,$(x),$*(x))], 

pn(v,/j,,T,<$>(x),<$>*(x)) = — r f dr [ dr\$(x)\ 2 - TrLnCT 1 + TrLnG 1 

V Jo J 



(2.7) 



(21 



is the one-loop effective action. This action includes in itself a series of terms containing 
derivatives with respect to $(x) and $*(x). In the lowest orders it corresponds to the 
Ginzburg-Landau effective action. 

The operation Tr in ( |2.8| ) is taken with respect to the space r, the imaginary time r 
and the Nambu indices. The action ( p.8|) is expressed through the fermion Green function 
which obeys the equation: 



t_$(t, r) + r + $*(r, r) 



G(t,t) = 6(t)6(t) 



(2.9) 
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with boundary condition 

G(t + /3,t) = -G(t,t). (2.10) 

The free Green function 

G (r,r) = G(r,r) 1*^=0 (2.11) 

in ( J2.8|) is needed to provide the regularization in the calculation of 
Q(v, T, <&*). The representation (|2.7|) , (|2.8|) is exact, although to perform the cal- 
culation in practice it is necessary to restrict ourselves to some approximation. Below 
we shall use the assumption, generally accepted in the 3D case (and in some quasi-2D 
systems), that the approximation including only the quadratic (Gauss) fluctuations of the 
fields <&(x) and <&*(x) (about their equilibrium values which will now also be denoted as 
$) describes the system quite well f\ 

The thermodynamical potential Q as well as the partition sum Z now depend on $ and 
$* which play the role of the order parameter. The order parameter appears due to the 
fact that one uses the effective potential Q(v , fi, T, $, <&*) instead of the exact potential 
Q(v, fi,T). Thus, it becomes necessary to write down the additional equations which 
determine the values of $ and $*. 

Moreover, in some cases it is sufficient to use the mean field approximation as is done 
in the original BCS theory. Here, however, we include the fluctuations also. Therefore, 
one obtains 



where 



Z(v, (i, T, $, <F) = exp[-f3n pot (v, /i, T, $, $*)] x 
J T?(A$)£>( A<F) exp j - J dr x dr 2 J dr 1 J dr 2 x 



/35 2 fi($,$*) 



5$*(ri,ri)5$(r2,r 2 ) 



A$(r 2 ,r 2 ) + 



<J> ,&*=const 



2 5§(t 1 ,t 1 )5$(t 2 ,t 2 ) 



A$(r 2 ,r 2 ) + 



<J> ,&*=const 



A$*(r 1>ri ) ~ 



/55 2 fi($,$*) 



2 5$*(ri,n)5$*(r 2 ,r 2 ) 



A$*(r 2 ,r 2 ) 



;2.i2) 



— [ dr\$(x) 



TTrLnG 7 + TTrLnG, 



-i 



;2.i3) 



is the mean field thermodynamical potential; A$(x) = $(x) — $ and A$*(x) = $*(x) — $* 
are the fluctuation deviations from the equilibrium value of the order parameter. 

In present review we restrict the consideration of the fluctuations to the critical line 
($ = $* = 0) only. So, the partition function ( |2.12p after integration over A$ and A$* 
acquires the form 



Z(v, fi, T) = exp[-f3Q pot (v, /i, T, 0, 0) - TrLnf 



-ii 



;2.i4) 



3 We note that this Gaussian approximation does not work in the 2D case at T ^ and, thus, we must 
modify it to analyze this case (see Chapter 5). 
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where now 



i-i, 



r, r = 



/M 2 fi($,$*) 



<5$*(r, r)5$(0,0) 



5(r)5(r) + tr [G(r, r)r_G(-r, -r)r+] |*,*. 



(2.15) 



is the inverse Green function of the order parameter fluctuations. 

To avoid possible misunderstanding, we shall write down the formulae for the Fourier 
transformations which are used throughout the paper; they connect the coordinate and 
momentum representations in the usual manner: 



F(iu n ,k) = J dr J drF(r, r) exp(iuj n T — ikr), 
+°° r (ik 

F(r,r)=T ^ / j—^F{iuj n , k) exp(-zo; n r + ikr), 

n=-oo J \ Zn ) 



(2.16) 



(2.17) 



where uj n = irT(2n + 1) are the fermion (odd) Matsubara frequencies and d is the dimen- 
sionality of the space (remember that in this Chapter d — 3). In the case of bosons the 
frequencies should be replaced by even ones: Q n = 2imT. 

For example, the Green function ( |2.9| ) has, in the momentum representation, the 
following form: 

iuoj + r 3 £(k) - r_$ - r+$* 



G(iu n , k) 



u 2 n +e(X) + \^\ 2 



(2.18) 



where £(k) = e(k) —jj, with e(k) = k 2 /2m and $ and $* are already taken to be constants. 



2.3 The mean field analysis 

Substituting ( |2.18| ) into fl2.13| ), one arrives at (see Appendix A) 
n pot (v,n,T, <!>,$*) = 

r i$i 2 



-2T 



dk 



2T 

where d = 3. 

The stationary condition 



V J (27r) d 

dk 



£ 2 (k) + |$| 2 
In cosh - — £(k) 



2T 



+ 



(27T) 



e(k) 

In cosh — — — e(k) 
2T v ' 



dn pot (v,^T c MF ^^ 



3$ 

results in the standard gap equation 

1 r dk 1 



V J (27r) 3 2£(k) 2T^ F 



tanh 



;2.i9) 



(2.20) 



:2.2V 
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Here we use the notation T c Mi? rather than T C 3D because Q pot describes the system in the 
mean field approximation only. In such a case T C MF does not depend significantly on the 
dimensionality of the space. Besides, as it will be seen, the value of T c Mi? may significantly 
exceed T C 3D (or the temperature of real condensation) in the strong coupling regime. 

Before proceeding further, we need to describe how we regulate the ultraviolet diver- 
gence in the gap equation ( 2.21|) . The idea is to replace the bare V by the low energy 
limit of the two body T-matrix (in the absence of a medium). In the 3D case we use the 



formula [[32 1 



rn 



1 

v + 



dk 



(o 22) 

47ia s V ' J (27r) 3 2e(k) 1 ' ' 

which defines the above-mentioned s-wave scattering length a s . As a function of the 
bare interaction, a" 1 increases monotonically from — oo for a very weak attraction to +oo 
for strong attractive interaction. Above the two-body bound state threshold in vacuum 
(aj 1 = 0), a s in fact is the "size" of this bound state with energy Ef, = —l/ma 2 s . The 
dimensionless coupling constant in the dilute gas model is then l//cpa s , which ranges from 
— oo in the weak coupling (BCS) limit to +oo in the strong coupling (Bose) one. 

Using ( |2.22|) and ( |2.21| ) one obtains directly the equation for the introduced transition 
temperature T C MF in terms of the renormalized coupling a s : 



m 



dk 



1 tanh^ 



2£(k) 



2T C MF 2e(k) 



(2.23) 



There are two unknown quantities in this equation T c and /i, and thus we need another 
equation, 



idn pot (v,^T c MF ^^*) 



n 



/• 



(2.24) 



which, as was already said, fixes the chemical potential for a given density. 

Note that unlike the BCS analysis in which one usually assumes that n — e F (the 
non- interacting Fermi energy), in the crossover problem fi will turn out to be a strongly 
dependent function of the coupling as one goes into the Bose regime where all the particles 
are affected by the attractive interaction. It directly follows from (|2.24|) and (|2.19|) that 
the number equation is given by 



MF\ 



n F (/i,T ( 



In the weak coupling limit, l/k F a s 
the well-known BCS result 



dk 



tanh 



71/ . 



(2.25) 



2T MF 

oo, we find from the system ( |2.23[ ) and ( |2.25| ) 



/i = e F , 



T, 



MF 



*0 



7T 



ire 



-e F exp 



2k F \a s 



;2.26) 



where 7 ~ 1.781. 

The equations can also be solved analytically in the strong coupling limit, where, 
however, one observes that the roles of the gap and number equations are reversed: the 
gap equation (|2.23j ) determines fi, while the number equation determines T C MF . In this 
limit l/k F a s — > +00 and one finds tightly bound (separate) pairs with the energy \E b \ — 
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l/ma 2 s (and \E^\ ^> ep). The non-degenerate Fermi system has here fi ~ — \Eb\/2 and 
its T C MF ~ \Eb\/2\n(\Eb\/eF) 3 / 2 . But really such a system can be hardly recognized as a 
Fermi system because all the fermions are bound in the pairs, so (at least at rather small 
temperatures) they in fact form Bose (local pair) system. 

This unbounded growth of the "transition temperature" is an artifact of the approxi- 
mation and there is, in fact, no sharp phase transition at T 6 MF (outside of weak coupling 
where T C MF = T c 3£) ignoring the small effects of thermal fluctuations). The point is that 
the mean field approximation becomes progressively worse with increasing coupling: the 
effective potential (|2.19| ) at $ = $* = can only describe a normal state consisting of 
essentially non-interacting fermions. 

While this is adequate for weak coupling, in the strong coupling limit unbound fermions 
exist in the normal state only at very high temperatures. In this limit, where the sys- 
tem is completely non-degenerate, a simple "chemical equilibrium" analysis (boson # 2 
fermions) yields a dissociation (pairing) temperature T^ ssoc = \Eb\/ ln(|i?&|/eir) 3 / 2 (~ Tp). 
We thus get convinced that for strong coupling T C MF is related to the pair dissociation 
scale rather than the T C 3Z) (-C T C MF ) at which the coherence is established. As it will be 
seen below the temperature Tp has an evident physical meaning in the 2D case: in the 
region Tbkt < T < Tp (further T p ) the system does not have a condensate (supercon- 
ducting state) but its one-particle spectrum acquires some features (in particular, the 
pseudogap) of superconductivity. 



2.4 Beyond the mean field approximation 



As we have seen the mean field thermodynamical potential (|2.13|) (or its direct form 
( |2.19| )) cannot be used to access the strong coupling limit. Therefore, to achieve this 
limit we have to include the effects of the order parameter fluctuations. It was already 
shown in Section 2.2 that treating the fluctuations at the Gaussian level is done by the 
partition function (|2.12|). It allows to incorporate the fluctuation into both the gap (|2.20|) 



and number (|2.24|) equations by adding the fluctuation correction to fl pot (v, /i, T, $, $*). 
As discussed in [|18l |32l in the 3D case it is sufficient, however, to take into account the 



fluctuation effects through the number equation only. That is why we do not need the 



correction d(£l — Vt pot )/d^\ 



to the gap equation. 



Consequently, we can restrict ourselves to calculation of the partition function at the 
critical temperature given by ( |2.21| ) that follows from (2.14) only. This permits us to get 



the fluctuation correction to the number equation (|2.24|) . Thus, one finally obtains from 



fl2T4) that 



n(v, ii, t, $ = $* = o) = n pot ( v , fi, t, o, o) + 

vT 00 r 
£ J dKlnr-^K) 



(2vr) ; 



(2.27) 



where 



r- 1 (zfi„,,K) = i 



dk 



2 J (27r) 3 £(k + K/2) +f(k-K/2) - iVt r 
" u £(k + K/2) £(k-K/2) 
tann h tanh ■ 



x 



2T 



2T 



(2.28) 
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Further, one has to remove the ultraviolet divergences from ( |2.28| ) by applying the same 
regularization procedure using the scattering length a s (see definition ( |2.22p ) as in case 
of the gap equation ( |2.21| ). With this substitution ( |2.28|) transforms to the form: 



i-i, 



zf2„,K) 



1 r dk 



3 



2 J (27,; 

unh £(k^ £(k-K^ 
2T 2T 



e(k) £(k + K/2) + £(k-K/2) -ifi, 

m 



■ x 



(2.29) 



According to [[TJ|, it is convenient to rewrite Q(v,fi,T, $ = $* = 0) in terms of a phase 
shift defined by T(u±i0, K) = \T(uj, K)| exp(±i5(uj, K)). Then the new number equation 
is given by the equality ("conservation law") 

n F (/i,T c 3D ) + 2n B (/i,T c 3D ) =n f , (2.30) 

with 

/ rr,3Ds f dK f°° du .d5{u,K) 

»-o.. T ' D ) S J^L aF"^-^ • (2 - 31) 

where ub(u)) = [exp(u;/T) — is the Bose distribution function. 

One sees from (|2.30|) that the system of fermions is separated into two coexisting 
and dynamically bounded subsystems: Fermi particles, or unbound fermions with density 
np(fi,T^ D ), and local pairs, or bosons with density n B (fi,T^ D ). Such an interpretation is 
possible and natural on the critical line only, when the finite density of condensate |$| is 
not included nonperturbatively in n^(/i, T). Besides, as we shall demonstrate in Chapter 
5, if the order parameter has a nonzero value, one can study the number equation without 
taking into account the bosonic contribution. 

The temperature Tj? D at which homogeneous long-range order is established is defined 
by the solution of ( |2.23p together with ( |2.3(J| ). At weak coupling (large nf), the results are 



essentially unaffected by the inclusion of Gaussian fluctuations in the number equation 
and T C 3Z) is the same as the mean field T C MF obtained above (see (|2.26|) ). 

In fact the equations can be solved in the Bose limit too. ^From ( |2.23| ) we find 
yu(T c 3D ) = —\Eb\/2, which is one-half the energy required to break a pair. Further, flOID 



can be simplified because the inequality \fi\/T 3> 1 is satisfied. From this inequality it 
follows that the isolated pole of T(u, K) on the real axis is situated far from the branch 
cut. The pole and the cut represent a two-body bound state with the center-of-mass 
momentum K and the continuum of two-particle fermionic excitations, respectively. The 
low energy physics for T T C MF is thus dominated by this pole, and one can approximate 
the phase at the pole by 5(u, K) ~ h6(uj — K 2 /4m + 2/i + \Ei\), so that 

nsQ*, Tf) = ^ / dKn B - 2„ - \E^ . (2.32) 



Finally, one gets 



/ \ 2 / 3 



T ' - {mm ™ = a218e - (2 - 33) 



which is simply the BEC result for bosons of mass 2m and density n//2. 
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Summarizing the underlying physics in these limiting cases it is necessary to stress the 
following [ IE ]: 

i) The value of T C 3D in the weak coupling limit results from thermal excitations of indi- 
vidual particles. 

ii) In strong coupling case, Tj} D results from thermal excitations of collective modes. These 
excitations are outside the range of mean field theory. That is one of the reasons why we 
took into account the fluctuations. 

The physics is thus quite different in the two limits considered: the pair breaking in the 
first case and the motion of bound pairs in the other. 

The results of the numerical solution of the equations are presented in [26j (see also 
p2| ). It may be worth commenting that the numerical result for T C 3D is a non-monotonic 
function of l//ci?a s with a maximum value at intermediate coupling which is slightly larger 
than the BEC value. The last is independent of the coupling constant. The situation is 
completely different for the discrete (lattice) model of the 3D crossover [IS] where Tj? D 
decreases as the coupling increases. So, certainly there is an optimal value of the coupling 
for T C 3D . The latter will also take place for the quasi-2D model which is in fact discrete 
in the third direction. We shall discuss such a model in Chapter 4. 



3 2D crossover: T = 

As we discussed in Introduction the 2D case has attracted much attention partly because 
its possible relevance to the layered HTSC. 

To avoid the above-mentioned problem of how to determine the T c value correctly in 
the strictly 2D system it is worthwhile to start with the zero temperature case. The prob- 
lem thus becomes, due to the integration over frequency (which replaces the summation 
over Matsubara frequencies), effectively a 3D one. Therefore, the 2D theorems M are not 
applicable to this case and one can speak about long-range order in the 2D systems at 
T = 0. 



We shall follow here the paper [24] (see also its extended version in [37| , pqj ) where the 
most general functional methods were used throughout. However, it is more convenient 
to go to the limit T = in the expressions from the previous Chapter rather than use the 



zero temperature technique as it was done in [24 . 



3.1 Model 



The model Hamiltonian that we shall consider coincides with ( [2 . 1| ) with one very im- 
portant exception: the dimensionality of the space here d = 2. This fact is crucial in 
the choice of the parameter which one has to change to trace the crossover. The fact is 
that 3D bound states in vacuum are known to form only if the corresponding coupling 
constant V exceeds some threshold. Thus, for the real cases one cannot achieve the Bose 
regime even at very low carrier densities if the attraction is not strong enough. This is the 
main reason why we have studied the 3D crossover as the function of the "renormalized" 
coupling a" 1 which corresponds to the "size" (radius) of the bound state at a^ 1 > only. 



A thresholdless bound state formation in the 2D space [39[ leads us to the important 
conclusion that one can reach the Bose regime by decreasing the density rif of bare 
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fermions at any coupling. So, to study the 2D crossover it is convenient and quite natural 
to regulate the density of carriers, or the Fermi energy e_p, which is the same for 2D metals 
with the simplest quadratic dispersion law: €f = nnf/m. In doing so, the coupling V 
should be replaced by its renormalized value which is the energy of the bound state in 
vacuum. We stress here once again that this parameter can be defined in the 2D case at 
any bare coupling V. The dimensionless coupling constant in the 2D case - the physical 
analog (see above) of l/kpa s in 3D - is thus given by the ratio e^/l^l, which changes 
from (in the Bose limit rif — > 0) to oo (in the BCS limit when rif — > oo). 



3.2 The Effective Action and Potential 



As we have already noted in the general case, (|2.8| ) is impossible for $ dependent on x. 
However, if one assumes the gradients of $ and $* to be small the action ( |2.8p can be 
naturally divided into kinetic and potential parts 

n(v,fi,$(x),$*(x)) = n kin (v,{i,$(x),$*(x)) + n pot (vifi, $,$*), (3.1) 

where the effective potential has been defined by ( |2.13| ) and calculated in Appendix A 
(its final expression is given by (|2.19|) ). The terms Qki n ($(x), $>*(x)) with derivatives in 
expansion ( |3.1p contains the important physical information, therefore we shall consider 
it in Section 3.4. 

Let us return to the effective potential. Going to the limit T — > one obtains from 

" <J>I 2 r Hk 

■ ' ' e 2 (k) + | $ |2_ e(k) 



V J (2tt) 2 

where the terms which do not depend on <3>, $* and fi are omitted. 

It is interesting that, in virtue of the invariance of the partition function 
respect to the phase transformation of the group 2/(1) 



(3.2) 



with 



-2ia 



$(x) 



$*(x) 



* f ( 



x)e 



$*(x)e 



2ia 



ID (see also (^TJ) and 



(3.3) 
can only be 



with real a the potential Q pot (v, $*) 
dependent on the invariant product $*$ 0. 

The analytic solution of the problem for the 2D case that we consider here is easier 
than the 3D one. Indeed, after performing the integration over k in ( |3.2| ) one directly 
obtains that 



+ 



W 



V An 



// 



In 



W-ij + J(W-ij) 2 + \<S>\ 



/i 2 + |$| 



w-fi + J(vy-/i) 2 + 1$| 2 v/x 2 + |$ 



(3.4) 



4 There is another transformation (when the si gn o f the phas e a is defined by the fermion spin rather 
than the charge) under which the Hamiltonian (2.1) (or (2.3)) is also invariant. Such a transforma- 
tion proves to be important for fermion- fermion repulsion (i.e. V < 0), or for the fermion-antifermion 
(electron-hole) channel of pairing. Apart from this difference the physics for the case of a repulsive 
interaction is identical to that under consideration. The complete set of gauge transformations for the 
Hamiltonian under consideration were originally given by Nambu [j3(| . 
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where the value W 
ary momentum. 



k. 2 B /2m is the conduction bandwidth and kg is the Brillouin bound- 



3.3 Main equations and analysis of solution 

If the quantity A is defined as the average value of |$|, then the equation for the 
extremum 

dVt pot {v, 



0. 



$=<J>*=A 



yields, according to 



A 



<9$ 



V 4^ 11 y/jj' 2 + A' 2 — jj 



(3.5) 



while the condition 



1 dtt P ot(v, /i, $,$*) 

V 



dfx 



n 



$=<J>*=A 

which sets the density of the particles in the system, takes the form 



W - J(W - fi) 2 + A 2 + Jfi 2 + A 2 = 2e F 



(3.6) 



(3.7) 



(3.f 



where we made use of the relation between e F and rif in the 2D case. 

Equations (|3.6|) and (|3.8|) , which were obtained in the mean field approximation which 
is quite sufficient at T = 0, make up a set for finding the quantities A and fi as functions 
of W and e F (or n/). It differs from the similar set of |2(J by the explicit dependence on 
W that, in principle, can be important for the case of narrow or multi-band systems P5" 



It should also be noted (see also the discussion after the equation (|2.24| )) that the 
necessity to utilize the system of equations in order to find self-consistently A and \i has 
been known for a long time (see |l|]). However, the fermion density in real 3D metals is 
virtually unchanged in practice so that, as a rule, the equation for \i is trivialized to the 
equation \x = ep and only the value A is regarded as unknown. The importance of the 
second equation for small particle densities was earlier pointed out in the papers |15|, [17| . 

Equations ( |3.6j ) and (|3.8|) allow, along with the trivial solution (A = 0, /i = ep) the 
nontrivial one: 



A 



e F (W-e F ) 
sinh 2 (27r/m\/) 



\i = e F coth 



2tt 



2tt 



W ( , 

coth ■ 

mV 2 \ mV 



(3.9) 



which is valid for any physically reasonable values of the relevant parameters. It is also 
very interesting that for small e F , there is a region where fi < 0, and that the sign 
changeover occurs at a definite point e F = W/2[l — tanh(27r/ml / )]. 

5 This is the parameter that is responsible for the appearance of a new (ordered, or with lowered 
symmetry) phase (see Section 2.2.). 
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The expressions that are found in [19|, [2(| follow directly from (|3.9| ) if, treating W as 



large and the attraction V as small, we introduce the 2D two-body binding energy 

e b = -2Wexp(-^j, (3.10) 

which does not include any many-particle effects. The introduction of the expression 
(|3.10D enables one to take the limit W 3> €f and thus justifies to a certain degree the 



use of the parabolic dispersion law. Second the fitting parameter e b is more physically 
relevant. For example it is well-defined even for potentials with repulsion. We stress that 
the introduction of e b instead of V also allows us to regulate the ultraviolet divergence, 
which is in fact present in the gap equation Q3.6|) . So, this step is quite similar to the step 
from equation ( |2.21 ) to ( |2.23| ). 



It should be mentioned that, in a dilute gas model, the existence of a two-body bound 
state in vacuum is a necessary (and sufficient) condition for a Cooper instability [19|, |2(| 



32]. This statement becomes nontrivial if one considers two-body potentials V(r) with 
short-range repulsion (e.g., hard-core plus attraction), so that one has to cross a finite (but 
really very weak f3~9f ) threshold in the attraction before a bound state forms in vacuum. 
Thus, by making use of ( |3.10|) it is easy to simplify (|3.9| ) |L9|, ^(J, p2| : 



A = y/2\e b \e F ; fl = - l -^ + e F , (3.11) 

with e F = \eb\/2. 

To understand the physical significance of these remarkably simple results we look at 
the two limits of this solution. For very weak attraction (or high density) the two-particle 
binding energy is extremely small, i.e. <C e F , and it is seen that we recover the 
well-known BCS results with strongly overlapping in r-space Cooper pairs. The chemical 
potential fi ~ e F , and the gap function A <C e F . 

In the opposite limit of very strong attraction (or a very low particle density) we 
have a deep two-body bound state \Eb\ 3> e F , and find that we are in a regime in which 
there is BEC of composite bosons, or "diatomic molecules". The chemical potential here 
fx ~ —\sb\/2, which is one half the energy of pair dissociation for tightly bound (local) 
pairs. 

It should also be kept in mind that in the local pair regime (/i < 0) the gap E gap 
in the quasiparticle excitation spectrum equals not A (as in the case fi > 0) but rather 
y/ fi 2 + A 2 (see [r?| and the review [p^l). 



Leaving aside the analysis of Vt pot for arbitrary values of the parameters (see [§7]), let 
us consider the most interesting case W — * oo, V — > with finite e^- Then finding from 
(|3.10|) the expression for An/mV and substituting it into (|3.4j) we obtain 



m 



fi 2 + |$| 2 - fi 




n pot (v,Li,$,<$>*) =«^|$r |^ v ' ' N ' : ===== ) • (3.12) 

whence we arrive at the expressions near point $, $* — > 0: 
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v^\ 2 In 



w 4¥ 



|$| 2 In 



l$l 



Equations (|3.6| ) and 



2^\e b \ 
can be written as 




(3.13) 



A 2 



^6 



/i 2 + A 2 + /i 



2e; 



(3.14) 



respectively; their solution, ( |3.11| ), is quoted above. As it is evident from ( p,13[ ), the 
potential term Q pot in the region /x < corresponds to particles with repulsion, which 
accounts for their BEC. We note here that for the 3D case this repulsion was obtained 
in [23 1 using a diagrammatic technique. In the region of high n/, where \x > 0, the Q pot 
cannot be represented as a series even for small $, which reflects the specific properties 
of the effective potential in 2D systems. 

If, proceeding from (|3.4j ) and ( |3.9|) , we find the thermodynamic potential difference of 
the trivial and non-trivial solutions, we can easily obtain the result that it takes form 



m 



tt pot (v, /i, A) - VL vot (v, fi, 0) = -v— 

In 



e 2 F coth 



2vr 
mV 



m 2 0(m) 



(3.15) 



which demonstrates (for various fi < < W/2) that at T = the nontrivial solution is 
always more favorable than the trivial one, and the point $ = $* = is unstable since 
here d 2 Q pot /d§d§* < for all allowed values of the parameters. 

As follows from (|3.12|) , both fl pot (/i, A) and Q pot (eF,0) are equal to —vmep/2ii, i.e. 
the potentials of the superconducting and normal phases turn out to be equal in the limit 
W — » oo, V — » 0. Nevertheless, it can be shown, just as in [p0| , that the superconducting 
phase has the lowest internal energy. In addition, the relevant difference is proportional 
to t>p(ep)A 2 (p(ep) = m/2Tr) being the density of states, which is energy independent in 
the 2D case); in other words, we arrive at the standard result of the BCS theory |], which 
is clearly valid for all values of the ratio e.p/|£&|. 



3.4 The gradient terms of the effective action 

Now we calculate the terms Qkin which contain the derivatives in expansion ( |3.1|) . As 
before, we shall assume the inhomogeneities of $ and $* to be small having restricted 
ourselves only to the terms with lowest derivatives. For simplicity we shall also consider 
the stationary case and calculate the terms with the second-order spatial derivatives only, 
which makes it possible to determine the coherence length £ and the penetration depth 
Xh of the magnetic field in a 2D superconductor. 

With these restrictions, taking into account the invariance of Q(v,fi, $, $*) (see equa- 
tion ( |2.8| )) with respect to the phase transformations ( |3.3| ) one can write a general form 
for the kinetic part of the action: 



Q Hn (v, fi, T, $(x), $*(x)) = T f dr J drT kin (<Z>, $*, V$, V<T) = 

(3.16) 



7' f dr I dr 



T 1 (|$| 2 )|V<I>| 2 + iT 2 (|$| 2 )(V|$| 2 ) 2 
2 
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where there are no items with a total derivative since the boundary effects are regarded 
as unessential, and the coefficients Ti i2 (|$| 2 ) are assumed to be unknown quantity. It 
follows from ( |3.16|) that in the second approximation in derivatives the variations in both 
the direction (phase) of the field $ and its absolute value are taken into account. 

To calculate the coefficients T% ^(l^l 2 ) we shall follow the paper |40| according to which 
one can determine the variation derivatives 



5 2 n km (<s> } & 



5$*(r)5$(0) 



5$(r)5$(0) 



$ ,&*=const 



<J>,<J>*=const 



-T f 'drlTAm 2 ) + |$| 2 T 2 (|$| 2 )]V 2 <5(i 
Jo 

-T f dT($*) 2 T 2 (m 2 W 2 5(r). 
Jo 



(3.17) 



Multiplying both equations (|3.17|) by r 2 and integrating over dr one can obtain 



drr 



2 5 2 n km (®^*) 



5$*(r)<5$(0) 



d ^ 2 5 2 Q km (^^) 



5$(r)5$(0) 



<£><!>* =CO?lSt 



<E><I>* =CO?lSt 



-AT / drtTxd^H + I^Tad^l 2 )]; 



AT / rfr($*) 2 r 2 (|$ 



(3.18) 



which allows one to determine the coefficients required. Indeed, let us define the correla- 
tion functions 



*y(q) 



dr exp(— iqr 



5 2 fW$,$*) 



«J$i(r)«J^(0) 



(3.19) 



$ ,&*=const 



It is then readily seen from (|3.18| ) and ( |3.19|) that Ti 2 (|^| 2 ) are determined from the 
derivatives of the correlators Zfy(q) in components of the vector q at q = 0: 



T 1 (|$| 2 ) = i 



T 2 (|$| 2 ) 



d 2 K 12 (q) $ d 2 K 22 (q) 



<9q 2 $ 
d 2 K 22 (q 



dq 2 



q=o 



4($ 



<9q 2 



(3.20) 



q=o 



it being apparently sufficient to know Ki 2 (q) and K 22 (q) only. Thus, the problem has 
been simply reduced to calculating the correlators. 

On the other hand, the effective action for the case of the imaginary time-independent 
fields can be written (see (|2.8|)) in the following form: 



n(v, //, T, $(r), $*(r)) = dr J dr {^|$(r) 



T ^ tr(r|Ln 



2m 



+ H] + r_$(r) + r + $*(r) 



|r) 



(3.21) 



where the operation tr refers to Pauli matrices, and the normalization term with Gq 1 is 
omitted. Comparing Q3.19Q and (|3.21|) and using the definition ( |2.9|) after the Fourier- 
transformation ( 2.17 ) we obtain directly: 

+0 ° r dk 

K 12 (q)=T 77rT^ tT [ G ( iuJ n,^)r + G(iuj n ,k + q)T_]; 

n=-oo J \ Zn ) 
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+°° r (fa 

K 22 (q)=T J j^ytr[G(iu} n ,k)T-G(iu; n ,k + q)' 



(3.22) 



where G(ico n , k) is the Green's function defined by (|2.18|) . Calculating the traces in (ft.22|) 
it is easy to arrive at the final expression for the correlators: 



nt^oo J ( 27r ) 

KM=T £ / 



K + ^( k )][o;2 + ^ (k + q)] 



(27T)2[^ + £ 2 (k)][^+£2 (k + q)] 



(3.23) 



where we used the notation £ 2 (k) = £ 2 (k) + |$| 2 . 

Summing over the Matsubara frequencies and then going to the limit T = one can 
obtain the final formulae: 



dk 



K 22 (q) 



2 £(k)+£(k + q) 
dk 1 



1 + 



k 2 

2^ ~^ 



(k + q) 2 u 



£(k)£(k + q) 



> ; 



*\2 



8tt 2 £(k) + £(k + q) S(k)£{k + q) ' 



(3.24) 



A somewhat tedious but otherwise straightforward calculation (see |37], |3~8"|) now yields 
for the desired coefficients of Q3.16Q : 



Ti(\®\ 2 ) 



1 fW-VL 



167T 



du 



2\§\ 2 -u 2 
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The expressions ( |3.25| ) and ( |3.26| ) together with ( |3.4| ) complete the calculation of all the 
terms of the effective action (R.ll) for the 2D metal with attraction between fermions. 



3.5 Correlation length and penetration depth versus doping 

Knowing T!(|$| 2 ) and T 2 (|<l ) | 2 ) one can find the values for different observables. For 
practical purposes it suffices to restrict ourselves to considering the coefficients obtained 
at the point |$| = A of the effective potential minimum. Besides, instead of ( |3.26 ) 
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it is convenient to introduce the combination T 2 (|$| 2 ) = T!(|$| 2 ) + 2A 2 T 2 (|$| 2 ) which 
determines the change in the |<3>| value only and arises as the coefficient at (V|$|) 2 . 

The analysis of (|3.25| ) and (|3.26|) reveals that at |$| = A both these functions are 
positive for ep < W/2 where they change their sign, which simply reflects the necessity 
to go over the antiparticle (from electron to hole or vice versa) picture for the region 
ep > W/2. On the other hand, in the region ep < W/2 the positiveness of these functions 
demonstrates the stability of the homogeneous ground state of the model concerned. 

One can succeed in simplifying (|3.25|) and ( p.26| ) if takes into account that, as a rule, 
e F < W (for example, in HTSC metal-oxides e F ~ 0). Then, using ( |Hq ) and 

(|3.11|) , one can rewrite Ti(A 2 ) and T 2 (A 2 ) in the very simple form: 
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where for completeness we have restored the Planck constant again. 

The explicit form of T\(A 2 ) and also a = «- 1 9 2 n pot ($, $*)/9$5$*||$| 2= 



(3.27) 



allows one 



to calculate the coherence length and the penetration depth. Since (see (|3.4j) ) 



m 



e F (W-e F ) coth^f 



* h W 2 coth 2 ^ 
mv 



in 



e F 



27rh 2 2e F + \e b \ 



- (W -2e F ) 2 
(W -> oo,V- 



0), 



(3.28) 



then, according to general theory of fluctuation phenomena | 41| , one gets that 
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This formula shows the dependence of £ on ep (or nj 
to compare (|3.29| ) with the definition of the pair size 
as the coherence length |H| |27|), namely 
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It is very interesting and useful 
p9|1 (also incorrectly referred to 



(3.30) 



where 



g(r) = ^(y BCS \4(r)4(0)\y B cs) 



(3.3i; 



is the pair-correlation-function for opposite spins and \^bcs) is the usual BCS trial func 
tion. For the 2D model under consideration the general expression (|3.30|) gives [20, 28 
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(3.32) 



where A and /x were given by ( p.ll|) . So, we are ready now to compare (|3.29| ) and ( p.32j ). 

Again to understand the underlying physics it is worth to look at two extremes of 
(|3.29|) and (|3.32|) . For high carrier densities ep 3> |e&| one finds that £ ~ £o ~ v F /A 



i.e. 
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the well-known Pippard's result is reproduced correctly. Moreover, if, according to fll9| , 
we introduce the pair size = h{m\eb\)~ 1 ^ 2 in vacuum, then it is clear that £ ~ £o^.£& 
when |£ft|<e_F. Therefore, both £ and £o prove to be close to the pair size which is much 
larger than the interparticle spacing. The latter statement one can see from the value of 
the dimensionless parameter t;k F ~ (,ok F ~ e F /A ;§> 1. 

In the opposite limit ^> e F of very low density one can see that <C £ ~ fc^ 1 
while £o ~ £&• Consequently the correct interpretation of £o is the pair size (in presence 
of the Fermi sea ( |3.30|) ) rather than the coherence length. The former in the extreme 
Bose regime is much smaller than the mean interparticle spacing, since £o&f *C 1, while 
always £k F ~ 1. The meaning of £ is indeed the coherence length because it remains 
finite and comparable with the mean interparticle spacing even when |e&| goes to infinity. 
This situation is consistent with the case of 4 He where the coherence length is nonzero 
and comparable with the mean inter- atomic distance although \e^\ (or energy of nucleon- 
nucleon binding) is really extremely large. Since ~ £o&f and since £bk F is directly 
related to the dimensionless ratio e F /\eb\ (which was discussed in Section 3.1) it can be 
inferred that ^kp is a physical parameter which can correctly determine the type of 
pairing. 

Calculations of the Xh value require one to introduce the external magnetic field H 
and also the usual extension of the derivatives TiW — > ftV — (2ie/c)A (e is the electron 
charge, c is the light velocity, A is the vector potential). Adding then the energy of the 
magnetic field 



dr I dr 



to the effective action 



Tj) by direct calculation one gets 
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Note that in ( |3.33| ) and ( |3.34j ) z is the number of superconducting layers per unit length 
taking into account the 3D character of the external field action. Hence, Ginzburg-Landau 
parameter (see ( |3.29| ), ( |3.34[ )) is equal to 
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Substituting here the accurate expressions ( |3.11fl , ( |3.27| ) and ( |3.28| ) for A, Tx(A 2 ) and a 
it is easy to find that ( |3.35| ) takes a simple enough form 
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(3.36) 



from which the explicit dependence of this parameter versus the ratio between e F and 
\eb\ follows immediately. So, for z ~ (10 7 — 10 8 )cm~ 1 (as occurs in real cuprates) the 
formula ( |3.36j ) for a large range of values for e F gives k ~ 10 2 . This means that 2D 
metals with low e F (so called underdoped case) turn out to be very strongly type-II 
superconductors. However, it follows from the same formula that with increasing e F the 
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parameter k decreases and in the limit e F ^> \e h \ (strongly overdoped case) it is possible 
(in principle) to obtain the value k<1 or in the other words, to change the type of the 
superconductor. 

We estimate finally the value of the second critical field H c2 . Indeed, equating the 
energy gain (see ( |3.15| )) of the state with A ^ and the field energy one can show that 



H, 



c'2 



For the same parameters we find that H C 2 ~ (10 5 — 10 6 )Oe = 10 — 10 2 T which also 
corresponds to the critical magnetic field values for the HTSC metal-oxides |42|, £43 



Finally, note that the concentration dependence of dH c2 /dT was studied also 44 



4 Crossover in the quasi-2D systems 

As we have already mentioned in the introduction there are many ways which allow one 
to extend the zero-temperature analysis of the 2D model from the previous Chapter to 
the case T ^ 0. Every way has some advantages as well as disadvantages. So, different 
generalizations of the model are to be investigated. Here we shall study a quasi-2D 
extension, which might be relevant for HTSCs with a relatively low anisotropy of the 
conductivity as, for example, is present in the Y-based cuprates. On the other hand, the 
quasi-2D model can be considered as an extension of the 3D model from Chapter 2. Our 



treatment is based on the paper |yj (see also its detailed versions [|45| ^q] ) 



4.1 Model 

The simplest model Hamiltonian density for the carriers in the quasi-2D system reads 



V 2 1 



_i_ 



cos (iV z d z ) + fi 



ip a (x) 



2m± m z d? z 

- V^(x)4(x)^(x)^(x), (4.1) 

where x = r, r±,r z (with r_i_ being a 2D vector); Tj) a (x), fi and V are already determined; 
m± is an carrier effective mass in the planes (for example, Cu02 ones); m z is an effective 
mass in the z-direction; d z is the interlayer distance. 

The Hamiltonian proposed proves to be very convenient for study of fluctuation stabi- 
lization by weak 3D one-particle inter-plane tunnelling. The extension of the 2D system 
in the third-direction could give rise to a possibility that it may undergo the 3D su- 
perconducting transition with the establishment of homogeneous long-range order. This 
transition does not contradict the 2D theorems [|J. We omitted in ([4.1|) the two-particle 
(Josephson) tunnelling considering it here to be smaller than the one-particle already 
taken into account. There can however be situations when Josephson tunnelling becomes 
more essential then the one-particle coherent one. In addition some authors consider that 
the most important mechanism for HTSC is the incoherent inter-plane hopping (through, 
for instance, the impurity (localized) states or due to the assistance of phonons.) 

It is significant that the large anisotropy of the conductivity cannot be identified with 
the similar anisotropy of the effective masses m z and m±. In particular, HTSC with 
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rather large anisotropy in z-direction do not display the usual metal behaviour at low 
temperatures. It means that interplane motion of the fermions is incoherent and the 
BKT transition could take place before the 3D superconducting transition (see the next 
Chapter). But, as will be seen, for justification of approximations used below the modest 
value of the ratio m z /m± > 10 2 is already sufficient. Such a value is present for instance 
in the HTSC YBa 2 Cu 3 6+ 5 (1-2-3). Because of this we shall study the simplest case of 
tunnelling which is described by the second term in ( |4.1| ). 

The method for the study of the Hamiltonian (|4.1|) mainly coincides with that used 
for the 3D model ( [2 . 1| ) . So, we shall often refer to the similar formulae from Chapter 2. 



4.2 The mean field approximation 



The effective potential ( |2.19| ) obtained in Appendix A is in fact model independent. Thus, 
one has to use it together with the dispersion law 



£(k)=e(k)-// 



2m 



m z d 2 z 



cos k z d z — ji 



(4.2) 



in accordance with the Hamiltonian ( |4.1| ). 

Again the stationary condition (|2.20 ) results in the gap equation ( 2.21) where now 



£(k) is determined by ( |4.2| ). 

It is necessary to emphasize that fot the temperatures of interest, the band width in 
the ^-direction is m~ 1 d~ 2 <C T C MF , i.e. the system is really a quasi-2D one. As for the last 
inequality, it is easy to see, that at m z « 10 2 m e and d z = 10A the value h 2 / \m z d? z kB) ~ 
10K is really far less then the usual critical temperatures in HTSC compounds. 

Now we need to describe how one can regulate the ultraviolet divergence in the gap 
equation ( |2.21| ). As has been already noted in Section 2.1 for the case of the local-pair 
superconductivity the cutoff of BCS type cannot be applied. Moreover, the regularization 
by the 3D scattering length a s (see the definition ( |2.22|) ) is unsuitable for the quasi-2D 
system. It turns out that in order to eliminate the divergences in ( |2.2ip it is convenient 
to introduce, as in the 2D case, the energy of the two-particle bound state (compare with 

e b = -2Wexp( ± , 4.3 

y m±V J 

where W = k. 2 LB /2m± is the bandwidth in the plane 0. The factor d z in the index of the 
exponent is necessary to preserve the dimensionlessness. 

Again, contrary to the usual BCS approach, the value of the chemical potential should 
be consistently defined from the equation ( p. 24 ), which leads to the second, or number, 
equation (|2.25| ). 

Taking into account the above-mentioned inequality m~ 1 d z 2 <C T c Mi? the number 
equation (|2.25|) in the limit W — > oo takes the following final form: 



2T MF ln 



1 + exp 



rpMF 



2e F , 



(4.4) 



6 Note that in the region of the parameters considered these bound states are formed without any 
threshold. 
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where (here) ep = 7mfd z /m± is the Fermi energy of free quasi-2D fermions with e(k) ~ 
k 2 - 

So, one needs to solve simultaneously the system of the equations ( |2.21|) with the 
dispersion law (|4.2j) and (|4.4j) with two unknown variables, T C MF and /i, respectively. 



At high carrier concentrations rif, such that ji ^> T C MF , the equality jj, ~ ep is indeed 
the solution of (|4.4j ). Than taking into account the regularization procedure, it follows 
from ( p.21| ) that (compare with ( |2.26| )) 

7 /^rn — 2 7 




This is just the well-known BCS result for a 2D metal 



In the opposite case of small n/, such as —(/, ^> T C MF , the roles of the gap ( |2.21| ) 



and number (44 ) equations are as above in some sense reversed: the equation fl2.21|) 



determines fi, while the (|4.4| ) determines the value of T C MF . Now, using the definition 
(|4.3Q one obtains from (|2.21| ) 

2W 1 r 27T r 2W 1 

In— = — / dt dx — . (4.6) 

\€b\ 2n Jo Jo x — 2ix + (2/m z dj ) cost 

Integrating the right side of (|4.6|) and in the limit W — > oo we immediately arrive at the 
final expression 

^ = ~~2~ " 2(m z diy\e b \- (4 ' 7) 

This expression, except for the second term, is identical to the result described in Sec- 
tion 2.3. The second term is directly connected to the quasi-2D character of the model 
and, despite its far smaller magnitude, is very important when the fluctuations are taken 



into account. At |ej| ^> ep the equation (4.4) transforms into the transcendental one: 



I I rpMF 

Ml ln^— . (4.8) 



2T C MF e F 

Its follows from (|4.8|) that we obtain a similar result to that in Section 2.3: for fixed 
ep, the value of T C MF grows rapidly as the coupling constant V increases. Thus, for the 
case of small carrier density, the temperature T C MF is not connected to the critical one 
T c , but in fact corresponds to the temperature of composite boson dissociation (compare 
with Section 2.3). 

4.3 The role of Gaussian fluctuations 

In order to investigate the effects of composite bosons formed in the system, one should 
again take into account the order parameter fluctuations. As was discussed in the previous 
Section, one can simply use the expressions from Chapter 2 with £ (k) determined by (|4.2j ) . 
Thus, the expressions ( 2.27 ) and ( |2.28| ) may be used directly if one applies the ultraviolet 



regularization by the energy (|4.3|) as has been done in Section 4.2. Once again, from 
( |2.27| ), we arrive at the "conservation law" ( |2.3(J| ). 
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Let us consider now the influence of n^^, T) determined by ( |2.31| ) on the behaviour 
of the system as the carrier density changes. We should solve self-consistently the system 
of equations (|2.21|) and ( |2.30| ). 

At high enough carrier density the contribution of ne(/i, T) to ( ^.3U| ) is negligible, and 
one arrives at the equality T c ~ T C A/F , where T c Mi? is given by ( f4.5| ). 

Recall (see Section 2.4) that in a more consistent scheme one should take into con- 
sideration the correction d(Q — Q pot )/d& to the gap equation ( |2.20| ), and consequently 
the fluctuations modify the equation. This correction may be especially important in the 
quasi-2D case. However, because of the conditions fi = ep ^> T c and m z d 2 z T c ^> 1, one 

fsee also 



can convince oneself that this correction changes T c rather weakly [[47 

Note that the consideration of the correction to the gap equation ( |2.20| ) should be 
very interesting because based on the general 2D theorems one knows that T c — > when 
m z — > oo (2D case) even though fi ~ ep. To trace this limit the above-mentioned 
correction must be taken into account. From the other side, when m z — > oo the BKT 
scenario will take place at any density rif and we have to modify the approach from 
Chapter 2 used here to study the formation of the inhomogeneous BKT phase. The latter 
will be accurately done in the next Chapter. 

Turning back to the equality T c ~ T C MF one can see, using ( |4.5| ) (the T C MF coincides 



with the expression obtained in |37] for a pure 2D superconductor due to the band nar- 
rowness in the /c 2 -direction), that T c increases with the growth of coupling constant V. 
This behaviour of T c is just as in the 3D case (see ( |2.26| )). 

At small concentrations, such that \fi\/T c ^> 1 taking into account the definition ( [4.3|) 
in the limit W — > oo one obtains from ( |2.28|) that 
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After introducing the phase 5(u), K) = — argT 1 (u + iO, K) it can again be written as: 



6(w,K) = n9[u 
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(4.10) 



At last, after substitution of the expression ( |4.7|) for the chemical potential into ( |4.10|) 
one arrives at the final equation for T c , namely: 



dK 



+ 
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cos K 7 d 7 
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(4.11) 



It is easy to see that the boson effective mass for its motion in the plane retains the 
value 2m^. As for the motion between the planes, the effective boson mass increases 
considerably: 2\eh\m 1 z d\{^> m z ). It can be stressed that this increase has a dynamical 
character as is testified simply by the presence of Physically, it is ensured by the 
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one-particle character (see equation (|4.1|) ) of the tunnelling between planes which only 
takes place through the virtual breakup of a pair for which the energy loss is of the order 

Now, using the formula for the BEC of an ideal quasi-2D Bose-gas [03 (see also the 



review |I5fl), it is easy to write down an equation for T c , which here takes the simple form 

T ~ T¥± = **. U 12) 

c 2m ± \n(2T c \e b \m 2 z di) 2 In (2T c \e b \m 2 z d*) ' v ' ' 

The last equation describes the characteristic properties of a quasi-2D superconductor 
with small carrier density: 

i) firstly, the critical temperature T c ~ ep, or (see above) T c ~ Uf/2, or the number of 
the composite bosons, as it should be in 2D case (recall that in a 3D one T c ~ n 2 / 3 (see 
(HH)), in contrast to the MF approximation 2D T C MF ~ ^nj > T c |g7| (see also (gill) 
because usually T C MF ~ A); 

ii) secondly, contrary to the case of the 3D superconductor where T c does not depend on 
V at all (see ( |2.33| )), in a quasi-2D system T c does depend on V, namely: T c decreases 



with the growth of V. As was stated above, the reason for this is the dynamical increase 
of the composite boson mass along the third direction. Thus, the growth of (or of 
V, which is the same) "makes" the system more and more two-dimensional even for the 
simplest case of a quasi-2D metal with a local four-fermion interaction. 

It is interesting to note that a decreasing T c can also take place in the case when 
the local pairs are bipolarons |50[| ; then, the increasing of coupling with phonons, which 



makes the pairs more massive, also leads to T c decreasing, (rather than increasing) as the 
electron-phonon coupling grows. 



The main results here are not only the expressions (O) and ( 4.12 ) for T c in the 
different limiting cases, namely the cases of large and small V. No less important is the 
comparison of these expressions which show that for a given density of fermions (i.e. a 
given ep) there are two essential regions. If \e b \ <C ep, then even in the case of small 
(by absolute value rif) densities the BCS formula is valid and T c grows with increasing 
\eb\ (see ( |4.5|) ). In the opposite case, \e b \ ^> ep, T c decreases with the growth of \e b \. It 
shows that in the case of quasi-2D systems (it seems that the HTSC belong to this case) 
T c (|£ft|) has a maximum. Consequently, there is a region (for fixed |e&|) of values of €f 
for which T c increases. If we also take into account the two-particle tunnelling in (|4.1| ). 
then the previous result will only be a lower bound for T c for large |ej,|. The region of 
e F ~ needs special study (see e.g. || 0]) because of the presence of both strongly 
developed fluctuations and the possible distinction of properties of such a Fermi-liquid 
from the Fermi-liquid of Landau type. 



5 2D crossover: finite temperature 

We have already studied the finite temperature (or more exactly T c ) crossover for the 3D 
and quasi-2D models. The 2D crossover at T = crossover was addressed also. So, we are 
ready to discuss the 2D crossover at T ^ 0. As was stressed in the Introduction (see also 
P2] |) an analytical treatment of the finite T crossover problem in 2D is quite a difficult task. 
This is primarily related to the necessity to treat T c as the BKT temperature T BKT below 
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which there is an algebraic (power decay) order and a finite superfluid density || (see 
also the review [55] and the book |53fl ). Nevertheless, some insight into the peculiarities 
of the formation of the BKT has been gained 0. 



5.1 Model and Formalism 



Most of previous analyses |54, 55, 56] of the behaviour of the 2D systems at T ^ have 



been based on a Nozieres-Schmitt-Rink approach As shown above, this is simply the 



Gaussian approximation to the functional integral which perhaps explains the difficulties 
faced in these calculations. On the one hand, Gaussian fluctuations destroy long-range 
order in 2D and if one searches for the T C 2D at which the order sets in one should get zero 
in accordance with the above-mentioned theorems ||. On the other hand, taking into 
account Gaussian fluctuations is completely inadequate to describe the BKT transition. 

Nonetheless, several steps have been made even in this direction. For example, the 
BKT transition in the relativistic 2 + 1-theory was studied in |57|] , and even the crossover 
from superconductivity to superfluidity was considered in |f22[ (see also ]|58[ ) according 



to the value of the carrier density nf (recall that rif = mep/n). However, the method 
employed in |22[ to the temperature Tbkt has a number of drawbacks. 



Specifically, the equation for T BKT was obtained neglecting the existence of a neutral 
(real) order parameter p, whose appearance at finite T, being due to the breaking of only 
a discrete symmetry, is consistent with Coleman-Mermin- Wagher-Hohenberg theorem J| . 
As we shall see below, p gives the modulus of a multivalued complex order parameter of a 
2D system as a whole, and only the modulus determines the possibility of the formation 
of nonuniform (including vortex) configurations in the system. 

However, as a result of allowing for a neutral order parameter, a region where p decays 
gradually to zero appears in the phase diagram of the system. This region separates the 
standard normal phase with p = from the BKT phase, where there is the power decay 
of correlations. Despite of the exponential decay of the correlations in it, this new region 
of states very likely possesses unusual properties, since p appears in all expressions in 
the same manner as does the energy gap A in the theory of ordinary superconductors, 
though to calculate the observed single-particle spectrum, of course, the carrier losses 
due to scattering of carriers by fluctuations of the phase of the order parameter and, in 



case of real systems, by dopants must be taken into account |59| . The possible existence 



of such a phase, which is also in some sense normal, might shed light of the frequently 



anomalous (see, for example, the reviews [34], |6(|) behaviour of the normal state of HTSC, 



specifically, the temperature dependencies of the spin susceptibility, resistivity, specific 
heat, photoemission spectra, and so on |H, 62], 21 1, for the explanation of which the idea 



of a pseudogap (and also spin gap) in the region T > T c is now widely employed. 

Thus, our objective in this Chapter is to calculate Tbkt and T p (T p is the temperature 
at which p — > 0) as functions of rif and to establish a form of this new region, which lies in 
the temperature interval Tbkt < T < T p . Besides, we will try to demonstrate using the 
example of the static spin susceptibility that this phase may really be used to explain the 
above-mentioned anomalous properties of the HTSC. That is the reason why the phase 
was called the "anomalous normal" phase. 

There is no need to write down the model Hamiltonian which is studied here, because 
it is identical with that of described in Chapter 3.1. 
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The desired phase diagram consisting from normal, anomalous and superconducting 
phases was calculated firstly in || employing the Hubbard-Stratonovich method (see 
Section 2.2, the equations (|2.6|) - (|2.8|) ). In the 2D case, however, instead of using the 
accepted method for calculation the partition function Z(v, /i, T, ) (see (PTB|) ) , one must to 
perform the calculation in modulus-phase variables. This prevents us from the subsequent 
treatment of the phase fluctuations at Gaussian level only. Thus, we will be able to take 
into account the phase degree of freedom with needed accuracy. 



The modulus-phase variables were introduced in accordance with ||63|| , where the 
parametrization 

$(x) = p(x)e- 2i0{x \ $*(x) =p(x)e 2W{x \ (5.1) 

was used. One can easily see that ( |5.1|) corresponds to the obvious transformation of the 
initial Fermi-fields, namely 

Mx) = X,{xY d{x \ V£(s) = Xl(x)e- i6{x \ (5.2) 

where the field operator Xv( x ) describes neutral fermions and exp[i9(x)] corresponds to 
the charge degree of freedom. In Nambu variables ( |2.2j ) the transformation (|5.2|) takes 
the following form 

= e iT3d{x) T(x) } ¥{x) = T\x)e- iTsd{x) . (5.3) 



Making corresponding substitutions ( |5.3j ) in the representation ( |2.6| ) and integrating over 
the fermi-fields T and we arrive at the expression (compare with ( |2.7|) and ( |2.8|) 



Z(v,fi,T) = J pVpVOexp [-/3n(v,fi,T,p(x),d9(x))\, (5.4) 

where 

PQ(v, fi, T, p(x), dO(x)) = ^ f dr [ drp 2 (x) - Tr LnG" 1 + Tr LnG 1 (5.5) 

V Jo J 



is as Q2.8Q the one-loop effective action, which, however, depends on the modulus-phase 
variables. The action ( |5.5| ) is expressed through the Green function of the initial (charged) 
fermions that has in the new variables the following operator form [] 



G- 1 = -Id T + r 3 + np(r, r) - 



V 2m J y 2m m J 

The free fermion Green function Go that provides the convenient regularization in the 
process of calculation was defined in (|2.11|) . It is important that neither the smallness nor 
slowness of the variation of the phase of the order parameter was assumed in obtaining 
expression (|5.6|). 



7 It may be obtai ned as a solution of some differential equation with the antiperiodic boundary condi- 
tions (see d]j|) and fl2.10p ). 
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Since the low-energy dynamics in the phases in which p ^ is determined mainly by 
the long- wavelength fluctuations of 6(x), only the lowest order derivatives of the phase 
need be retained in the expansion of Q(v, p, T, p(x), d9(x)): 



Q(v, p, p(x), dO{x)) ~ fL k in(v, p, T, p, dd(x)) + Q pot (v, p, T, p) 



where 



n km (v,p,T,p,d6(x)) =TTrJ2 -(G^) n 



n=l 



n 



p=const 



and (see (|JJ)) 



n 



pot 



[v, p, T, P)=(h J drp 2 - TTr LngT 1 + TTr LnG 1 



p=const 



(5.7) 
(5.8) 

(5.9) 



The kinetic flkin and potential Vt pot parts are expressed through the Green function of the 
neutral fermions which obeys the equation (compare with (|2.9|)) 



g{r,r) = 6(r)8{r) 



(5.10) 



and the operator 



E(d6) = r 3 id T 9 



2m 



2m 



m 



(5.11) 



The representation (|5.7| ) allows one to get a full set of the equations which are necessary 
to find out T BKT , p(T BKT ) and p(T BKT ) at given ep (or, for example, p(T) and p(T) at 
given T and €f). While the equation for Tbkt will be written using the kinetic part ( |5.8| ) 
of the effective action, the equations for p(Tbkt) and p(Tbkt) (or p(T) and p{T)) could 
be obtained using the mean field potential ( |5.9| ). It turns out that in the phase where 
p 7^ the mean-field approximation for the modulus variable describes the system quite 
well. This is mainly related with a nonperturbative character of the Hubbard- Stratonovich 
method, i.e. most of effects are taken over by a nonzero value of p. 



5.2 The equation for T, 



BKT 



If our model were reduced to a some known model describing the BKT phase transition, 
we would be easily write the equation for T BKT . Indeed, in the lowest orders the kinetic 



term ( |5.8| ) coincides with so called classical XY model |52|, |53| which has the following 
continuum Hamiltonian 

H = J - Jdr(V6(r)) 2 . (5.12) 

Here J is the some constant (in the original classical discrete XY model it is the value of 
spin) and 9 is the angle (phase) of the two-component vector in plane. 

The temperature of the BKT transition is, in fact, known for this model, namely 

Tbkt = \j- (5.13) 
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Despite the very simple form of the equation (|5.13| ), it was derived || (see also |52 



|5^|) using the renormalization group technique, which takes into account the non-single- 
valuedness of the phase 9. Thus, the fluctuations of the phase are taken into account at 
a higher approximation than the Gaussian one. 

To expand Qkin up to ~ (V#) 2 , it would be sufficient to restrict ourselves to terms with 
n = 1, 2 in the expansion (|5.8| ). The procedure of calculation (see Appendix B) is similar 



to that employed in pl |, where the case of large densities rif at T = was considered, 
and gives f\ 



n 



T rP 



kin 



2 Jo 



dr I dr 



where 



J(/i, T, p(/i, T))(V6>) 2 + K(ji, T, p(/i, T))(d T 6) 2 
J(fi, T, p) = —n F (p, T, p) / dx j 

IJl 7T J-U/2T „ / 



cosh 2 W x 2 + ^ 



(5.14) 



(5.15) 



characterizes the stiffness of neutral condensate 

m i _ p 



K(/,,T,p) = -ll + 7 = f = f tanh 



IT 



(5.16) 



and a value 

n F (p,T, p) 




p 2 + p 2 + p + 2T\n 



1 + exp 



T 



(5.17) 



has a sense of fermi-quasiparticles density (for p = the expression ( |5.17| ) is simply the 
density of the free fermions). Note that J(p, T, p = 0) = 0. 

A direct comparison of the expression ( |5.14|) with the Hamiltonian of XY model ( |5.12 ) 
makes it possible to write a equation for Tbkt- 
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J{h,Tbkt, p(^,Tbkt)) — Tbkt- 



(5.18) 



Although mathematically the problem reduces to a well-known problem, the analogy is 
incomplete. Indeed, in the standard XY model (as well as the nonlinear a model) the 
vector (spin) subject to ordering is assumed to be a unit vector with no dependence 
on T. In our case this is fundamentally not the case, and a self-consistent calculation of 
Tbkt as a function of nf requires additional equations for p and p, which together with 
(|5.18|) form a complete system. 

8 The exponentially small correction is omitted here. 
9 A total derivative with respect to r is omitted. 

10 There is no doubts that in certain situations (for example, very high T) it also can become a ther- 
modynamical variable, i.e. dependent on T, as happens in problems of phase transitions between ordered 
(magnetic) and disordered (paramagnetic) phases when the spin itself vanishes. Specifically, for quasi-2D 
spin systems it is virtually obvious that as one proceeds from high-T regions, at first a spin modulus 
forms in 2D clasters of finite size and only then does global 3D ordering occur. 
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5.3 The effective potential and the equations for p and /i 



There is no need to repeat the calculation of the effective potential. The point is that 
the effective potential ( |2.19| ) calculated in Appendix A depends on (see Section 3.2) the 
invariant product $$* = p 2 only. Thus, one may immediately write 



tt P ot(v,H,T, p) 



El 
v 



dk 



2T In cosh 



e(k)+p 2 



2T 



-£(k) 



(5.19) 



where £(k) = k 2 /2m— p. Then the desired missing equations are the condition dVL pot (p)/dp 
that the potential ( |5.19| ) be minimum and the equality v^dVtpot/dp = —rif, which fixes 
rif. For them we have, respectively: 



1 

V 



dk 



(27r) 2 2 ^2( k ) +p 2 



tanh 



e(k)+p 2 

2T 



(5.20) 



n F (p,T,p) = n f . (5.21) 

The equations (|5.20|) and ( 5.21|) obtained above comprise a self-consistent system for 
determining the modulus p of the order parameter and the chemical potential p in the 
mean-field approximation for fixed T and rif. 

As we have already discussed in Section 3.3, the energy of two-particle bound states 
6b (see its definition ( |3.10| )), is more convenient to use than the four-fermi constant V. 
For example, using the identity 



i x 
tanh — 

2 



1 - 



one may easily go to the limits W 
this renormalization becomes 



exp x + 1 

oo and V — > in the equation ( |5.20| ), which after 



In 



\ £ b\ 



\/p 2 + p 2 - p 



du 



-H/T 



u 2 + 



'R)2 



exp yV + (^) 2 + 1 



(5.22) 



Thus, in practice, we will solve numerically the system of the equations ( |5.18| ), ( f5.22j ) and 
(|5.21|) to study Tbkt as function of rif. 

It is easy to show that at T = the system ( |5.22| ), ( |5.21| ) transforms into the system 
(|3.14|) which was already addressed. Recall that its solution is p = J^e^e? and p = 
— \eb\/2 + €F- This will be useful for studying the concentration dependencies of 2A/Tbkt 
and 2A/T p , where A is the zero-temperature gap in the quasiparticle excitation spectrum 

ill 

[ VP? + p 2 , p<0. K ' 

Setting p = in the equations (|5.20 ) and ( ^.21| ), we arrive (in the same approximation) 
at the equations for the critical temperature T p and the corresponding value of p: 



ln N7 

T p 71 



v/ 2T p tanhw 

du 

o u 



(7 = 1-781), 



(5.24) 
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T„ln 



1 + exp (¥. 



(5.25) 



Note that these equations coincide with the system which determines T^ 2D ^ MF and n(Tj; 2D ^ MF ) 
(see ||37|| ). This is evidently related with mean- field approximation used here. There is, 
however, a crucial difference between these values. Namely, if one takes into account the 
fluctuations, the value of T C 2D should go to zero, while the value of T p should stay finite. 
That is the reason why the temperature T p has its own physical meaning: the incoherent 
(local or Cooper) pairs begin to be formed formed below T p . At higher temperatures there 
are these pair fluctuations only (see e.g. 



5.4 The phase diagram 



The numerical investigation of the systems ( |5.18| ), ( |5.22| ), ( |5.21| ) and ( |5.24j ), ( |5.25| ) gives 
the following very interesting results, which are displayed graphically. 
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Figure 1: Tbkt and T p versus the free fermion density. The dots represent the function 
p(€f) at T = Tbkt- The regions of the normal phase (NP), anomalous normal phase 
(ANP) and BKT phase are indicated. 

a) The anomalous phase region (see Fig.l) in the present model is commensurate with 
the BKT region. But it has not been ruled out that in the case of the quasi-2D model this 
region will disappear as rif increases. For example, in the case of an indirect interaction 
it was shown |65| that the anomalous phase region really exists at the low carrier density 
only, i.e. it shrinks when the doping increases. 

b) For low €f(<^. \sb\) the function Tbkt(^f) is linear, as it also confirmed by the analyti- 
cal solution of the system ( |5.18p , ( |5.22| ) and ( |5.21| ), which gives Tbkt = ep/2. Surprisingly 
such behaviour (although for T c ) is copied in different families of non-conventional ("ex- 
otic") superconductors (HTSC including) with comparatively small Fermi energy [[2~l|j . 

We note that in this limit the temperature T c of formation of a homogeneous order 
parameter for the quasi-2D model (see Chapter 4, the equation ( 4.12j )) can be written in 
the following form |65| (see also f2T|| ) 



J-c ~ 



T 



BKT 



2\ ' 



\n(e F \e b \/4tf,) 



(5.26) 
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where t\\ = l/{m z d 2 z ) is the inter-plane hopping (coherent tunnelling) constant. This 
shows that when T c < Tbkt the weak three-dimensionalization can preserve (in any case, 
for low n/) the regions of the anomalous and BKT phases, which, for example, happens 
in the relativistic quasi-2D model \ ffl>\ . At the same time, as the three-dimensionalization 
parameter tu increases, when T c > Tbkt the BKT phase can vanish, provide, however, 
that the anomalous phase region and both temperatures T p and T c are preserved. It 

follows from ( 5.26Q that the BKT phase vanishes when t\\ > J2\£{,\eF{= A). 
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Figure 2: fi(T) for different values of e F /\e b \: 1 — 0.05; 2 — 0.2; 3 — 0.45; 4 — 0.6; 5 - 
1; 6 — 2; 7 — 5. (For jj, > and \x < the chemical potential was scaled to e F and |£&|, 
respectively.) The thick lines bound regions of the BKT, anomalous normal and normal 
phases. 

c) Figure 2 shows the values rif for which \i differs substantially from e F and, in other 
words, the Landau Fermi-liquid theory becomes inapplicable for metals with low or inter- 
mediate carrier density. As expected, the kink // at T — T p , experiments on observation 



which were discussed in |B7J and have been interpreted for the 1-2-3 cuprate |5S| , becomes 
increasingly less pronounced as increases. But in the present case it is interesting that 
in the approximation employed it happens at the normal-anomalous phases boundary or 
before superconductivity really appears. Therefore it would be of great interest to perform 
experiments which would reveal the temperature dependence /x(T) especially for strongly 
anisotropic and relatively weakly doped cuprates. 

d) It follows from curve 3 in Fig. 2 that the transition (change in sign of fi) from local to 
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Cooper pairs is possible not only as ep increases, which is more or less obvious, but also 
(for some nj) as T increases. 




2 3 



Figure 3: 2A/Tbkt and 2A/T P versus the free fermion density. 



however, somewhat lower and reaches the 
Eb\. It is interesting that this concentration 



e) Finally the calculations showed (see Fig. 3) that the ratio 2A/Tbkt is always greater 
than 4.4. The value 2A/T p (= 2A/T C MF ) is 
BCS theory limit of 3.52 only for ep ^> 

behaviour is consistent with numerous measurements of this ratio in HTSC 
that the divergencies of 2A/Tbkt and 2A/T P at ep — > are directly related with the 
definition (|5.23|) . 
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Discussing the phase diagram obtained it can be emphasized that the qualitative 
ideas about the crucial importance of phase fluctuations for underdoped HTSC was, to 
our knowledge, first discussed by Emery and Kivelson ||. These authors, starting from 
the experimental data of Uemura et al. 



21 and some others (see the Refs. cited in 



and from the well-known general observation the superconducting state is characterized 
by a complex order parameter, introduced the temperature Tg at which phase order can 
occur. They have also argued that at low superconducting carrier density and poor 
screening (i.e. in bad metals) phase fluctuations become more significant than all other 
fluctuations, so that the classical XY-model is suitable for underdoped cuprate oxides. 
When T c <C T™ ax (T™ ax is the temperature at which the phase order disappears if the 
disordering effects of all the other degrees of freedom are ignored), phase fluctuations 
must be relatively unimportant and the observable T c will be close to the T C MF predicted 
by BCS-Bogolyubov theory. Otherwise, T c ~ T^ ax < T C MF , and T C MF is simply the 
temperature at which local pairing occurs. 

Undoubtedly, the approach developed in this Chapter is self-consistent and a more- 
or-less complete (in the hydrodynamical approximation) extension of the semiqualitative 
results presented in |], At the same time there is only one difference (but in our 
opinion an essential one) between the behaviour of the function we obtain for T c MF (ny) = 
T p (rif) in the limit rif — > and the function sketched for T c Mi? in ||. It is clear from the 
figure (see Fig.l) (and can be shown analytically) that in our case the zero density limit 
is T p (0) = 0, in contrast to Refs. || ^TJ where T C MF grows when rif decreases. 

On the other hand, the limit T p (rif — > 0) — > (the same argument also applies to 
Tdissoc in Section 2.3) cannot be considered sufficiently regular due to the strong increase, 
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for small rif, of the neutral order parameter fluctuations which were not taken into account 
in the approximation used. From the physical point of view perhaps the most consistent 
limit for this, extremely low (when /i < 0) fermion density, region is: T p (jif — > 0) — > Tp « 
2 v /p T T7? ~ \e b \ (see ( ^23; ) and Fig. 2). 



5.5 "Spin-gap" behaviour in the anomalous normal phase 

It would be very interesting to study how a nonzero value of the neutral order parameter 
affects the observable properties of the 2D system. Does this really resemble the gap 
opening in the traditional superconductors, except that it happens in the normal phase? 
Or, in other words, does the pseudogap open? 

We shall demonstrate this phenomenon taking, as a case in point, the paramagnetic 
susceptibility of the system 

To study the system in the magnetic field H one has to add the paramagnetic term 



H 



PM 



p B H \ip\(r)^(r) - V>{(r)V>i(r) 



(5.27) 



to the Hamiltonian ( |2.1| ). Here p B = eh/2mc is the Bohr magneton. Note that, using 
the isotropy in the problem, we chose the direction of field H to be perpendicular to the 
plane containing the vectors r. (Recall that in this Chapter d = 2 and r is the 2D vector.) 
It is a very simple matter to rewrite TipM m the Nambu variables ( |2.2|) , namely 



n 



PM 



-IM B H& '(r)tf(r). 



(5.28) 



Then, adding the corresponding term to the equation (|5.10|) for the neutral fermion Green 
function, it is easy to obtain that of in the momentum representation (compare with(|R2|)) 



Q(iu n ,k,H) 



(iu n + p B H)I - r 3 £(k) + np 
(iu n + p B H)I + r 3 £(k) - Tip 



(tu n +p B H)2-e(k)-p 2 ' 

The static paramagnetic susceptibility is expressed through the magnetization 

dM(p,T,p,H) 



(5.29) 



x(t*,T, p) 



OH 



(5.30) 



H=0 



which in the mean-field approximation may be derived from the effective potential 

w/ rp m l dn pot (v,p,T, p,H) 
M(p,T,p,H) = -- — . 



(5.31) 



Thus from ( 5.31]) one obtains Q 

oo 

M(p,T,p,H)=p B T £ 



n=— oo 



dk 



tr[g{iuj n ,k,H)I}. 



(5.32) 



11 That was done with V.P. Gusynin. 

12 Note that one should use the first part of flS.29 ) to get this formula. 
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Then using the definition ( |5.30| ) one arrives at 

dk 



x(j*,T,p) = p B 



2T y: 



e 2 (k)+ P 2 -^ 



The sum in ( |5.33|) is calculated in Appendix C (see the equation 
obtain the final result 



1 r 

X(^,T,p) = XPauli^ J_ 



dx 



M// cosh 2 Jx 2 + -§f2 



(5.33) 
)) and, thus, we 

(5.34) 



where XPauii = f i2 B m / 7T is the Pauli paramagnetic susceptibility for the 2D system. 

To study x as a function of T and rif (or ep) the formula ( |5.34|) should be used together 
with the equations Q5.22 ) and (|5.21|) . 

For the case of the normal phase (p = 0) one can investigate the system analytically. 



Thus ( |5.34 ) takes form 



x(^,T,p = 0) = xp. 



1 



auli 1 / i rn \ j 

1 + exp[—p/l ) 

while p is already determined by (|5.25 ). This system has the solution 

x(e F , T, p = 0) = XPauit[l - exp(-e F /T)], 
which coincides with that known from the literature 0. 
1 
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Figure 4: X {T) for different values of e F /\e b \: 1 — 0.2; 2 — 0.6; 3 — 1; 4 — 5; 5 — 10; 6 
- 20. 



The results of numerical study of the system (|5.34j) , ( |5.22|) and ( |5.21| ) are presented in 
Fig. 4. One can see that the kink in x happens at T = T p as in the dependence of p on 
T. Below T p the value of x(T) decreases, although the system is still normal. This means 
that the spin-gap (pseudogap) opens. The size of the pseudogap region depends strongly 
on the doping (e.p/|£&|), as it takes place for the real HTSC |JT], [6(J. For small values 
of ej?/|£b| this region is large (T p > 2T BKT ), while for the large ratio it is small. 
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6 Concluding remarks 



To summarize we have discussed the crossover in the superconducting transition between 
BCS- and Bose-like behaviour for the simplest 3D, quasi-2D and 2D models with s-wave 
direct nonretarded attractive interaction. It has been pointed out that optimally doped 
HTSC are still on the BCS side of this crossover, although they are certainly far away 
from the standard BCS description. 

Above we emphasized the model description which in our opinion proves to be the 
most suitable for the clarification of the diverse physical properties peculiar to electronic 
systems with changeable carrier density in any dimension. While there is still no gener- 
ally accepted microscopic theory of HTSC compounds and their basic features (including 
the pairing mechanism), it seems to us that this approach, although in a sense phe- 
nomenological, is of great interest since it is able to cover the whole region of carrier 
concentrations (and consequently the whole range of coupling constants), temperatures 
and crystal anisotropies. It, as we tried to demonstrate, allows one not only to propose a 
reasonable interpretation for the observed phenomena caused by doping but also to predict 
new phenomena (for example, pseudogap phase formation as a new thermo dynamically 
equilibrium normal state of low dimensional conducting electronic systems). 

Evidently there are a number of important open questions. They may be divided into 
two classes: the first one concerns the problem of a better and more complete treatment 
of the models themselves. The second class is related to the problem of to what extent 
these model are applicable to HTSC compounds and what are the necessary ingredients 
for a more realistic description. 

Regarding our treatment of the 3D and quasi-2D models, it is obvious that one has 
to take into account the interaction between the bosons also. In particular the Gaussian 
approximation is not sufficient to give reliable results for T c at intermediate coupling 
when the "size" of the bosons is comparable with the mean distance between them. As 
for the 2D model, there is some unconfirmed numerical result J7TJ] based on a fully self- 
consistent determination of a phase transition to a superconducting state in a conserving 
approximation, which state that the superconducting transition is not the BKT transition. 
Besides, it would be very interesting to obtain the spectrum of the anomalous normal 
phase. 

Concerning the question to which the models considered are really applicable to HTSC, 
it is obvious that most of the complexity of these systems is neglected here. We did not 
take into account an indirect nature of the interaction between the fermions and d-wave 
pairing. Note, however, that some attempts to study the crossover for these cases were 
made [72], [73], ^8] (for general review see |3JJ]). 

A lot of peculiarities of HTSC are now connected with stripe structure of Cu02 planes, 
which, according to many experiments (see [74]] and references therein), are divided onto 



bands of the normal and superconducting bands. There is interesting and important 
problem how to investigate these systems. 

The problem of the crossover from BEC to BCS (especially for 2D systems) is such 
a rich one that without doubt will bring us a lot of surprises in the near and far future. 
One of them is perhaps the unified theory of superconductivity and magnetism [75] which 
has already excited a lot of interest and criticism (e.g. [76|). 
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A The effective potential 



Let us derive the effective potential (|2.f9|) . To obtain it one should write down the formal 
expression Q2.13| ) in the momentum representation, so that 



Q pot {v,fi,T,$,<f>*) =v 



V 



+00 



t E 



dk 

W) 1 



+ T E 



dk 



■tr[lnG l (iuj n , k)e 



■tr[lnG 1 {iuj n ,k)e 



, 5^+0, 



i5u> n T 3 ] 



where 



and 



G L (iu n , k) = iujJ - r 3 £(k) + r+$ + r_$* 



G^^Umk) = G 1 (iu n ,k) 



' iu n - f (k) $ 

$* iu n + f (k) 

iu n — e(k) 



(A.l) 



(A.2) 



(A.3) 



$=$*=^=o I £<j n + e(k) y 
are the inverse Green functions. The exponential factor e lSuJnT3 is added into (|A.1|) to 



provide a right regularization which is necessary to perform the calculation with the 
Green functions (see 17711). For instance, one obtains that 



lim^ E tr[lnG (iu n ,k)e 



iSuj n T 3 ] 



n=— 00 
+00 



,5^ m 1 E t r [ m C 1 (^n, k)] cos(W n 4- 



^ n=— 00 



i E sin 5d; n tr[(ln G *(zo; n , k) — lnG {—iuj n: k))r^[ 

+ OO 



E trpnCT^k)]- 



£(k) 



T 



where the properties that 



InG (iu n ,k) = -t 3 



e(k) 



and 



were used. 



E 



sin 8u r . 



OJn>0 Un 



' sindx 1 . . 

ax = -— signo 

2vrT Jo x AT 6 



(A.4) 



To calculate the sum in (|A.4j) one have to firstly use the identity tr In A = In det A, so 
that ( |A.1| ) takes the following form 



n pot (v,fi,T,$,<5>* 



V 



+00 



- t e 



<ik j detG 1 (zo; n ,k) 



cik 



(2tt) c 



-£(k) + e(k)] 



(A.5) 
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Calculating the determinants of the Green functions (|A.2 ) and ( |A.3| ) one gets 



n pat (v,tM,T,<f>,$*) 



^ + e 2 (k) + |$| 



(A.6) 



where the role of Go(iu n , k) in the regularization of Q pot is evident now. The summation 
in (|A.6|) can be done if one uses the following representation 



In 



u 2 n + a 2 



dx 



u)\ + a 2 + x cu 2 .+b 2 + x / 



Then, the sum 



E 



i 



, :0 (2k + l) 2 + c 2 4c 
may be now applied and one obtains 



7T TIC 

tanh — 
2 



In 



cu 2 + 6 2 



1 



tanh 



2y/WTx 2T 2v / ^ 2 + x 

Integrating ( A.9 ) over x one thus arrives to the expression: 



tanh 



2T 



r E 



dk ln ^ 2 + e(k) 2 + i$i 2 



(2vr 



^ 2 +^ 2 (k) 



2T 



dk , cosh[^ 2 (k) + |$| 2 /2T] 
In - 



{2n) d cosh[e(k)/2T] 
Finally, substituting flOCD into (TOQ we get ( g]Tg ). 



(A.7) 



(A.8) 



(A.9) 



(A.10) 



B Low energy kinetic part of the effective action 

Here we derive the kinetic part ( p.8| ) of the effective action ( |5.7| ). To obtain it one should 
calculate directly the first two terms of the series in ( |5.8|) which are formally written as 

2- 

rP r T ~ r ( (W) 2 



n2. = T J dT J dr J^~2 Y,J dktr[g(iu; n ,k)T 3 } (i^ + L-i-j , (B.l) 



f2fc2i = TTr(£/£) and fij^ = |TTr((?£(?£). The straightforward calculation of fljj.^, gives 

dktr[£(iu; n , k)r 3 ] ^ic 

/*/ • i \ + T 3^(k) - Tip 

g{tu n ,k) = — t^t— 2 B.2 

is the Green function of the neutral fermions in the frequency-momentum representation 
(compare with (|2.18|) ). The summation over Matsubara frequencies u n = n(2n + 1)T and 



where 
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integration over k in ( p.l|) can be easily performed using the sum ( |A.8| ) and thus one 
obtains 



n { £ n = T J drj drn F (ji,T,p) lid T 9 + 



2m 



(B-3) 



where n F (p,T, p) is determined by ( 5.17 ). We note that S has the following structure 
S = r 3 Oi + /O2 where 0\ and O2 are some differential operators (see ( |5.11|) ). One can 
see, however, that the part of E, proportional to the unit matrix /, does not contribute 

For the case T = |63|, |64| when real time t replaces imaginary time r, one can argue 
from the Galilean invariance that the coefficient of dt9 is rigidly related to the coefficient 
at (V#) 2 . So it does not appear in f^in- We wish, however, to stress that these arguments 
can not be used to exclude the appearance of the term (V6>) 2 from Vt^ n when T^O, thus 
we must calculate it explicitly. 

The 0\ term in £ yields 



T rP 



dr I dr 



T 



(2vr) 2 n 



00 . 

/ d\ttx[G(iLO n i\t)TzG{iu ni \t)T2\ X 



id T 9 + 



2m 



(B.4) 



And from ( |B.4|) we find that 

^l(Oi) = ~ [dr J dvK(p,T,p) (id T 6 + 



{V0f 
2m 



(B.5) 



where K(p,T,p) was defined in ( |5.16|) . It is evident that 0\ term does not affect the 
coefficient of (V#) 2 . Further, it is easy to make sure that the cross term from 0\ and O2 
in &fa n is absent. Finally, the calculations of the 2 term contribution to £l\ in Q give 



/ / dr I dr 



T 



(2tt) s 



00 „ 

J dkkhr [g (iu n , k) ig (iu n , k) i] x 



(W 

4m 2 



;b.6) 



Thus, after summation over Matsubara frequencies (see Appendix C) 



dr / dr- 



1 



32ir 2 m 2 



dk- 



cosh 



e(k) + p 2 



,(vef 



(B.7) 



2T 



As expected this term vanishes when T — > but at finite T it is comparable with 
Combining ([BJ), (|B~5|) and we obtain (p34]) . 
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The higher than (V6 1 ) 2 derivatives were not found here. 
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C Summation over Matsubara frequencies 

Here we perform the summation over Matsubara frequencies in the following expression 

oo 

T Yl tr[G(iw n ,k)ig(iu n ,k)t\, (C.l) 

n=— oo 

where the Green function Q{iuj n ,\s) is given by (B.2j ). At first, using the elementary 
properties of the Pauli matrices one obtains 

tr[£(^„, k)ig(iu nj k) J] = g^M + l^ . (C.2) 



Then, the summation can be easily carried out, if one uses the following sums WE 



V" ^ n na n 2 1 

^ [(2fc + l) 2 + a 2 ] 2 = 8^ ^ T " 16a 2 cosh 2 ^ ( } 

and 

~ (2fc + l) 2 _ 7T TTa 7T 2 1 

^ [(2k + l) 2 + a 2 ] 2 ~ 8a t& 2 + 16 cosh 2 ™ ' 1 ' 

Assuming that a 2 = (£ 2 (k) + p 2 )/n 2 T 2 one directly arrives to the final result 

~ £ 2 (k)+p 2 -^ 2 1 1 fc 

n £^ 2 + £ 2 (k)+P 2 ] 2 ~ 2T 2V / e(k)+p2 ' l ' j 

COsh ^rp 
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